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Martensites are long-lived nonequilibrium structures pro- 
duced following a quench across a solid state structural tran- 
sition. In a recent paper (Phys. Rev. Lett. 78, 2168 (1997)), 
we had described a mode-coupling theory for the morphol- 
ogy and nucleation kinetics of the equilibrium ferrite phase 
and twinned martensites. Here we calculate nucleation rates 
within a first-passage time formalism, and derive the time- 
temperature-transformation (TTT) diagram of the ferrite- 
martensite system for athermal and isothermal martensites. 
Empirically obtained TTT curves are extensively used by 
metallurgists to design heat treatment cycles in real mate- 
rials. 

PACS: 81.30.Kf, 81.30.-t, 64.70.Kb, 64.60.Qb, 63.75.-fz 

A martensitic transformation |^ is a nonequilib- 
rium solid state structural transition which rcsuhs in a 
metastable phase {martensite) , often consisting of aher- 
nating twinned arrays within the parent sohd (austenite). 
A characteristic of this transformation is that it is diffu- 
sionless, i.e., the velocity of the transformation front is 
much larger than typical diffusional speeds of atoms. In 
contrast, infinitesimally slow cooling results in the equi- 
librium phase (ferrite). 

Martensites are commercially important, for example 
as the hard constituent in steel, or in shape-memory al- 
loys like Nitinol. These desirable properties, depend on 
the amount of transformed martensite. Over the years 
a tremendous amount of empirical and theoretical 
[H,D knowledge regarding the kinetics of this transfor- 
mation has accrued. A convenient representation, used 
extensively by metallurgists is the time-temperature- 
transformation diagram Q]. The TTT diagram (Fig. 
1), is a family of curves parametrized by the fraction 
6 of transformed product. Each curve is a plot of the 
time required to obtain 5 versus temperature of quench. 
It is clear from the figure that the shape of the curves 
are qualitatively different for the ferrite and martensite 
products. There are two important features regarding 
the TTT curves for martensites. («) The transformation 
curves lie below a temperature M^, the martensite-start 
temperature. This means that a martensite does not ob- 
tain at temperatures above Ms, even if the quench rates 
are infinitely high, (ii) The curves are parallel to the 
time axis. This implies that the martensitic transforma- 
tion is completed immediately following the quench (such 
martensites are called af/ierma^ martensites). Such trans- 
formation curves have been obtained empirically ; in the 



absence of a unified nonequilibrium theory for marten- 
site and ferrite nucleation, these curves have not been 
'derived'. 

In a recent paper j|] we presented a mode-coupling the- 
ory for the nucleation and growth of a product crystalline 
droplet within a parent crystal. We showed that for 
slow quenches, the droplet grows diffusively as an equilib- 
rium ferrite inclusion, while for fast quenches, the droplet 
grows ballistically, as a martensite having twinned inter- 
nal substructure, with a speed comparable to the sound 
velocity. Given this unified description, can we arrive at 
the general features of the phenomenologically obtained 
TTT curves ? The details of the underlying physics are 
ofcourse complex — this includes an understanding of 
the heterogeneous nucleation and growth of martensitic 
grains and other intervening phases, and the subsequent 
collision of grains emanating from correlated nucleation 
events. In this Letter, we demonstrate that certain qual- 
itative apsects of the TTT diagram can nevertheless be 
obtained using very general features of the distribution 
of nucleation events and the nucleation dynamics arising 
from our mode-coupling theory ^ . 
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Fig. 1 Experimental TTT curves (adapted from Ref. Q]) 
for the structural steel AISI 1090, containing approximately 
0.84% C and 0.60% Mn. The letters A, F and C represent 
the f.c.c.-austenite, the equilibrium b.c.c-ferrite and the car- 
bide precipitate (FeaC) respectively. Curves corresponding 
to 0, 50 and 100 % transformation are shown. Below a tem- 
perature Ms, the metastable martensite (M) is formed - the 
transformation curves for martensites are horizontal. 
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Though our analysis can be generahzed to any struc- 
tural transition in arbitrary dimensions, we will, for cal- 
culational simplicity, focus on the square to rhombus 
structural transition in 2-dimensions. It will be clear 
that the qualitative form of the leading 'edge' of the TTT 
curves (labelled 0% in Fig. 1) is independent of dimen- 
sionality and the type of structural transition. 

The square to rhombus structural transformation in- 
volves a shear -|-volume deformation, and so the strain 
order parameter has only one nontrivial component 
e — {dyUj. -\- dxUy)/2 (u is the displacement vector field). 
A quench across this transition results in the nucleation 
of a rhombic product in a square parent. At the ini- 
tial time, the transformed region is simply obtained as 
a geometrical deformation of the parent creating a local 
atomic mismatch This leads to a discontinuity in 
the normal component of the displacement field across 
the parent-product interface Au • n (n is the unit normal 
to the interface) This discontinuity appears as the 
vacancy field £,4> = £.in-int — n^acj/n, where ^ is the in- 
terfacial thickness, n is the average number density, and 
riint and riyac are the interstitial and vacancy densities 
respectively (measure of the compression or dilation of 
the local atomic environment) ||^. 

The free-energy functional J- describing this inhonioge- 
neous configuration contains the usual bulk elastic free- 
energy of a solid Fei and an extra interfacial term de- 
scribing the parent-product interface In our 2-d ex- 
ample, the bulk free energy F^i is constructed to have 
three minima — one corresponding to the undeformed 
square cell (e = 0) and the other two corresponding to 
the two variants of the rhombic cell (e = ieg). The total 
(dimensionless) free-energy functional to leading order in 
6 is, 
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{Vef + l{^dr,efe- (1) 



The modulus 7 = VIq~^Y^' ^'^dndnc{r) (prime denotes a 
sum across the interface) is the surface compressibility of 
the vacancy field and depends on the local orientation of 
the parent-product interface (c(r) is the direct correlation 
function of the liquid at freezing). The three minima 
of the homogeneous part of F at e = (square) and 
e = ±eo = ±[(1-|-Vl — 3a)/3]"'"/^ (rhombus), are obtained 
in the parameter range < a < 1/3. The first-order 
structural transition from the square to rhombus occurs 
at a = 1/4. For 1/4 > a > 0, the square is metastable. 
The degree of undercoohng 1/4— acxTj — T, where Ts 
is the temperature at which the equilibrium structural 
transition occurs. 

A quench across the structural transition, nucleates a 
region of the product (e = ±eo) within the parent (e = 0). 
The growth of this nucleus is described by a Langevin 
equation for the broken symmetry variable e and the va- 
cancy field (p. Moving with the growing interface, the 
equations of motion become, 

■V (2) 



90 „ „9<5J^ 



(3) 



where the noise ry is an uncorrelated gaussian 
white noise with variance proportional to temperature 
{r]{x,t)r]{x',t')) = 2kBTr6{x-x')S{t-t'). is the mi- 
croscopic vacancy diffusion coefficient which has an Ar- 
rhenius dependence D0 = Z3oo exp(— A/ZcsT). Typical 
values of parameters for pure (99.98 %) Fe at around 
1223 - 1473 °K have been measured to be H A = 
280 - 310 kJ-mol-i and = 0.4 - 4.0 cmVs, which 
fixes ~ 10~^^ cm^/s. 

Our free-energy functional naturally admits two widely 
separated time scales — r;, the relaxation time of the 
order parameter to the local minima and t„, the first- 
passage time for the order parameter to go from the 
local to the global minimum. The shorter time scale 
T; ^ [^"(e = 0)]~^^^ — ^/c (c is the velocity of trans- 
verse sound) with a typical value of lO^^'^s The re- 
laxation time of (f>, given by ~ ^^/D^, lies between 
Ti <^ Tfj, < Tn- Thus the order parameter is 'slaved' to 0, 
which acts as a time dependent source coupling to 9„e. 
The 'instantaneous' value of 4> is given by the Eq. (|^). 

The exact formula for the nucleation rate of a solid 
nucleus growing within a solid parent is unknown. In a 
continuum field theory approximation, the first passage 
time r„ would be given by the Kramers formula with 
a time-dependent energy barrier j9j. In this paper we 
shall take t„ = F"^ exp(Ai?*), where /S.E* is the energy 
barrier (in units of ksT) to form the critical nucleus. We 
have ignored the curvature dependent prefactors since 
they are only weakly dependent on (j){t) and a. 

To calculate the energy of the critical nucleus, we use a 
variational ansatz for the strain profile defining a nucleus 
of linear dimension i, 



Co if -L/2 + ^/2<x<L/2-^/2 
and -L/2 + ^/2<y < L/2 - ^/2 

if -L/2-^/2>x> L/2 + ^/2 
or -L/2-^/2>y> L/2 + £,/2 



(4) 



Within the interface of width ^, the strain e linearly 
interpolates between and eo (see Fig. 2(a)). The va- 
cancy field (j){x,y,t) is obtained from a solution of the 
diffusion equation Eq. (||) , with the initial condition that 
4>{x, y, 0) = Au • n, where Au can be computed from the 
variational ansatz Eq. (m ; 



[x, y, t) = $(x, y - L/2, t) + $(y, x - L/2, t) + 
$(x, y + L/2, t) + $(y, x + L/2, t) 



(5) 



within the interface and without, where 



dt 



L/2 .C/2 

$(a,b,t) = "1^ / dx' / dy'x'x 
{a-x'f + {b^y'f\ 
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(6) 
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Fig. 2 (a) Variational ansatz for e showing a cross -section 
across a square shaped nucleus of size L and interface width ^. 
(b)-(d) 3-d plots showing a sequence of 4>{x,y,t) surrounding 
the nucleus at different times. 

Figures 2(b)-(d) is a sequence of three-dimensional 
plots of (t>{x,y,t) at different times. A positive (nega- 
tive) <f) indicates a higher interstitial (vacancy) density. 
The free-energy of a nucleus of size L at a given time 
t, E{L,t]S^), is obtained by substituting the variational 
form for e and the 'field' (j){x,y,t) in Eq. (|l|). For every 
L, t, we minimize the energy E(L, t ; ^) with respect to the 
variational parameter ^, to obtain E(L,t). The energy 
AE* and size L* of the critical nucleus for every time t, is 
determined by dE{L,t)/dL = 0. The energy of the 

critical nucleus decreases with time, and so nucleation of 
the stable phase will occur when the time dependent nu- 
cleation barrier becomes small enough. This is obtained 
by self consistently solving 



T„ = r-i cxp(A^*(r„)). 



(7) 



Once the nucleus exceeds the critical size it grows. 
Thus the leading edge of the TTT curves, defined as the 
time required to form (5-^-0+ amount of the product at 
a given temperature of quench, should be given by Eq. 
(0). This however does not take into account the rate 
at which nuclei are produced. There are principally two 
kinds of nucleation modes [|lO| — (i) Athermal and (ii) 
Isothermal. A number of alloys like the one portrayed in 
Fig. 1, follow the athermal mode of martensite produc- 
tion. The amount of martensite formed is a function only 
of the temperature of quench and not the time of holding 
at that temperature pT] . A convenient parametrisation 
is 



where to Tn and No{T) increases with decreasing T be- 
low Ms- It is clear that as long as to ^ Tn, all grains are 
essentially nucleated at once, and the leading edge of the 
TTT curves is given by r„ (solid line in Figure 3). The 
curves in Figure 3, have been derived using the quoted 
values for the parameters Doc and A for Fe alloys [ p2[ . 
The modulus 7 is fixed by the condition that the inter- 
cept CM, (7) = Ms/iTs. For Fe alloys, Ms = 100-200°C 
and Ts = 900°C, leading to a 7 = 0.1. The time axis is 
in units of F^^ taken to be 10~^s. All these parameters 
can be varied over a wide range without changing the 
qualitative nature of the results. 

For low undercooling (small 1/4— a), the critical barrier 
height is large. The cf) field relaxes quickly, and so r„ 
asymptotes the — 00 curve (dotted line in Fig. 3). 
Once the critical nucleus is formed, it grows as a ferrite 
For larger undercooling, decreases reaching a 
vanishingly small value for a < aM^- This part of T„(a) 
asymptotes the = curve, when the (f) field remains 
frozen at the parent/product interface. As shown in Ref. 
Q , the nucleus twins in the direction of motion resulting 
in a martensite. 
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Fig. 3 Calculated leading edge of the TTT curves showing 
ferrite and athermal martensite regions (solid line). The val- 
ues for Doo and A in dimensionless units have been translated 
from the quoted values for Fe (see text) — Doo — 4x 10^'* and 
A — 8. The arrow indicates the martensite-start temperature 
hi our units. The dashed curve shows the first passage 
time for = 00 

Our calculation thus successfully reproduces the two 
distinctive features of athermal martensites — the hori- 
zontal martensite transformation curves, and a well de- 
fined martensite-start temperature which is indepen- 
dent of Doo and decreases with increasing 7. 

On the other hand in isothermal martensites, the nuclei 
are thermally generated, and so to is large compared to 
Tn- In this situation, the leading edge of the TTT curve 
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does not have a precise meaning. We can nevertheless de- 
fine the leading edge from the product (dN/dt) . Be- 
ing produced by a thermal activation process, dN/dt de- 
creases with increased undercooling. However as stated 
above, the critical barrier height and hence t„ decreases 
with increased undercooling. This produces the charater- 
istic C-shaped TTT curves of the isothermal martensite 
10 . The qualitative features of the leading edge do not 
change if the isothermal nucleation is autocatalytic. 

The other curves in the TTT family (Fig. 1) corre- 
spond to the transformation of a fixed amount of the 
product. This requires a detailed knowledge of the dy- 
namics of patterning generated by the nucleation and 
growth of grains , together with the shape and size of 
the individual grains. 

In summary, we have derived the essential features of 
the leading edge of the time-transformation-temperature 
curves for both athermal and isothermal martensites, 
within the first-passage time formalism. This makes use 
of a recently developed [|j mode-coupling theory for the 
morphology and nucleation kinetics of the equilibrium 
ferrite and twinned martensites. Such transformation 
curves have thus far been obtained empirically, and to 
the best of our knowledge have never been derived. 

Most martensites contain alloying elements which oc- 
cur as interstitial or substitutional impurities. The sol- 
ubility of these impurities may be different in the par- 
ent and product phases (eg., interstitial carbon in steel). 
This would induce large scale diffusion of impurities as 
the transformation proceeds (eg., carbon diffuses away 
from the ferritic product) . TTT diagrams of such alloy 
steels often show a secondary bulge together with the for- 
mation of an intermediate structure called bainite. These 
issues may be understood within our formalism by cou- 
pling the density of impurities to the strain e. We inves- 
tigate the role of alloying in a forthcoming publication. 

We thank K. P. N. Murthy for a discussion on time- 
dependent nucleation barriers. 
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